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The energy dependence of multiplicity fluctuations was studied for the most central Pb + Pb 
collisions at 20A, 30A, 40A, 80A and 158A GeV by the NA49 experiment at the CERN SPS. The 
multiplicity distribution for negatively and positively charged hadrons is significantly narrower 
than Poisson one for all energies. No significant structure in energy dependence of the scaled 
variance of multiplicity fluctuations is observed. The measured scaled variance is lower than the 
one predicted by the grand-canonical formulation of the hadron-resonance gas model. The results 
for scaled variance are in approximate agreement with the string-hadronic model UrQMD. 
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1. Introduction 

At high energy densities (w 1 GeV /fm 3 ) a phase transition from hadron gas to quark-gluon- 
plasma (QGP) is expected to occur. There are indications that at RHIC and top SPS energies quark- 
gluon-plasma is created at the early stage of heavy ion collisions |[l], g]. The energy dependence 
of various observables show anomalies at low SPS energies which might be related to the onset of 
deconfmement [^j). Lattice QCD calculations suggest furthermore the existence of a critical point 
in the phase diagram of strongly interacting matter which separates the line of the first order phase 
transition from a crossover. Models predict an increase of multiplicity fluctuations near the onset 
of deconfmement [Q] or the critical point 

This motivates a vigorous experimental and theoretical study of multiplicity fluctuations in 
high energy nuclear collisions In addition to the effects mentioned above also the "back- 

ground" statistical fluctuations are interesting to study. In a grand-canonical statistical model scaled 
variance co l is close to one. In the past it was commonly believed that the result of the grand- 
canonical, canonical and micro-canonical ensemble should be similar in the infinite volume limit, 
like it is for mean quantities. Recent calculations showed that this assumption is wrong. The scaled 
variance in the infinite volume limit is largest for the grand-canonical ensemble and smallest for 
the micro-canonical one [||, ^]. It is because the fluctuations in the canonical and micro-canonical 
ensembles are reduced by conservation laws and this reduction is, unlike for the mean, volume 
independent for V — > °°. 

The first preliminary results on the energy dependence of multiplicity fluctuations in central 
Pb + Pb collisions at SPS energies were presented in this paper in order to look for experimen- 
tal signatures of increased fluctuations due to a phase transition or the critical point and for the 
predicted reduction of fluctuations in relativistic hadron gas due to conservation laws. 



2. The NA49 Experiment 

The NA49 detector [10] (see figure [|) is a large acceptance fixed target hadron spectrometer. 
Its main devices are four large time projection chambers (TPCs). Two of them, called vertex TPCs, 
are located in two superconducting dipole magnets with a total bending power of 7.8 Tm. The 
other two TPCs are installed behind the magnets left and right of the beam line allowing precise 
particle tracking in the high density region of heavy ion collisions. The measurement of the energy 
loss dE jdx allows particle identification and a good rejection of electrons in a large momentum 
regime. 

A large variety of different colliding systems was studied. Beams of p and Pb are available at 
the CERN SPS directly, C and Si beams were produced via fragmentation of the primary Pb beam. 
Targets of liquid hydrogen or solid foils of different materials were used. Data on central Pb + Pb 
collisions at 20A, 30A, 40A, 80A and 158A GeV were recorded and analysed for this analysis. 



'a commonly used measure of fluctuations, CO = " , where Var(n) and < n > are variance and mean of multi- 
plicity distributions, respectively. 
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Figure 1: Setup of the NA49 experiment and different target configurations, a) for Pb + Pb, b) for p+p and 
c) for p+A collisions. 
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Figure 2: Collimator in front of the Veto calorimeter at 158A GeV [pT|. 



2.1 Centrality Determination 

The downstream veto calorimeter, originally designed for NA5, allows a determination of the 
centrality of a collision by measuring the energy in the projectile spectator region [jljj]. A collimator 
is located 25 m downstream from the target and is adjusted for each energy in such a way that all 
projectile spectator protons, neutrons and fragments can reach the veto calorimeter. For 158A GeV 
the hole in the collimator is with respect to the beam axis ±5 cm in vertical direction and —5 cm 
and +38 cm in horizontal direction taking into account the deflection of charged particles by the 
magnetic field (figure Q). For 40A and 80A GeV the hole is ±12 cm in vertical and —13 cm 
respectively +47 cm in horizontal direction. For 20A and 30A GeV the ring calorimeter with a 
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Figure 3: Acceptance of the Veto calorimeter for neutral (left) and positively charged (right) main vertex 
particles at 158A GeV as a function of total momentum p and transverse momentum pj. 

hole of radius 28 cm and positioned 18 m away from the target and 10 cm away from beam axis in 
horizontal direction is used as a collimator. For higher energies the ring calorimeter is positioned 
17 cm off in horizontal direction and has only a small influence on the acceptance of the veto 
calorimeter. 

For the following analysis the collisions are selected according to their energy in the veto 
calorimeter. A large fraction of the Veto energy is due to projectile spectators. However also a 
fraction of non-spectator particles, mostly protons and neutrons, contribute to Veto energy. The 
acceptance of the Veto calorimeter at 158A GeV is shown in figure || for neutral and positively 
charged particles. 

The centrality C of an event with a veto energy Ey eto can be calculated using the known trigger 
centrality C m - g = ^ L and the veto energy distribution as: 

C _ °Eve, _ r ^ V "" dN I dE Vet0:tri g 

— ^trig ' poo ... I j j-, v^- 1 / 
Ginel Jo "A / aEyetofrig 

where dN /dEv e to,trig is the Veto energy distribution for a given trigger. 



2.2 Track Selection 

Since detector effects like track reconstruction efficiency might have a large influence on mul- 
tiplicity fluctuations, it is important to select a very clean track sample for the analysis (figure |]). 
Therefore the acceptance for this analysis is limited to a part of the forward hemisphere, where the 
NA49 detector has the highest tracking efficiency. This was done by restricting the analysis to the 
rapidity interval 1 < y(n) < y heam 2 for 20A to 80A GeV and 1.08 < y{n) < 2.57 for 158A GeV 3 . In 



addition a cut on transverse momentum according to [ ]13fl was applied, which is dependent both on 
rapidity and azimuthal angle. This ensures that only tracks with sufficient points both in the vertex 
TPC 2 and in one of the main TPCs are used. For them the reconstruction efficiency is larger than 
98%. The disadvantage of this track selection is that the fraction of accepted tracks is small and 
depends on collision energy, see figure B. 



2 Rapidity is calculated in the center of mass system assuming pion mass. 

3 A different rapidity cut is used for this energy because of missing geometrical acceptance near beam rapidity [jl4j]. 
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Figure 4: Acceptance for negative hadrons at 158A GeV as a function of rapidity (pion mass assumed) and 
Pt (left), as well as a function of (j) and pj (1.4 <y < 1.6, right). 
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Figure 5: Acceptance for negative hadrons as a function of collision energy calculated using Venus model 
and GEANT detector simulation. The acceptance was also estimated using the ratio of measured multiplicity 
of negatively charged hadrons to the published or preliminary NA49 total n~ and K yields, the difference 
(< 5%) is shown by vertical error bars. 



A cut on the distance between the interaction point and the track trajectory extrapolated to the 
target is used to reduce the background of secondary particles originating from weak decays (track 
impact parameter) 4 . Removing this cut would change the scaled variance by less than 2%. 
A cut on the energy loss of particles in the detector gas (dE /dx) removes electrons and allows the 
study of hadron multiplicity fluctuations only. The scaled variance will be less than 2% different if 
this cut is removed. 

3. Fluctuation Analysis 

The multiplicity distribution in central (7 — 10%) Pb + Pb collisions at 158A GeV is shown in 
4 The difference has to be smaller than 4 cm orthogonal and 2 cm parallel to the magnetic field direction. 
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figure |6|. The measured distribution is significantly narrower than a Poisson one. 

The basic measure of multiplicity fluctuations used in this analysis is the scaled variance: 

co = ^l (3.1) 

< n > 

where Var(n) and < n > are variance and mean of multiplicity distributions, respectively. 

The scaled variance for positively (co(h + )), negatively (oo(h~)) and all charged hadrons (co(h ± )) 
will be presented. The scaled variance of a Poisson distribution is 1, independent of its mean mul- 
tiplicity. A larger CO might indicate additional non-statistical fluctuations, a smaller CO might be a 
hint for a suppression of fluctuations e.g. due to conservation laws. 

Figure ^ shows the centrality dependence of multiplicity fluctuations in Pb + Pb collisions 
at different energies for central collisions. In general CO decreases with increasing centrality (i.e. 
decreasing Veto energy), this trend is stronger for higher energies. 

The results for collisions selected by minimum bias and central on-line trigger are shown 
separately in figure [7] and they agree within errors. For further analysis the results for central data 
sets are used because of larger statistics in comparison to the minimum bias data sets. 



3.1 Experimental Biases 



In superposition models (e.g. the Wounded Nucleon Model [12]) the total number of particles 
produced in a collision n is the sum of the number of particles produced by different independent 
particle production sources: 

n = Y J n i,so (3.2) 

i 

In these models the scaled variance has two contributions. The first is due to the fluctuations of 
the number of particles emitted by a single source co so , the second is due to the fluctuations in the 
number of sources C0u '. 

co = co so + < n so > co k (3.3) 
where < n so > is the mean multiplicity of hadrons from a single source. 

The number of sources related to projectile participants is fixed by fixing the energy detected in the 



Veto Calorimeter and consequently the number of projectile spectators (see chapter |2.1[ ). 

For very central collisions analysed here the Veto energy is dominated by the energy of forward 
going participants and not spectators. For these collisions the number of projectile participants is 
in good approximation fixed and independent of Veto energy. Therefore no correction for the finite 
width of the Veto energy bins is applied. 

The finite resolution of the Veto calorimeter introduces additional fluctuations in Np r0J even for 
a fixed Veto energy. For very central collisions their influence on the scaled variance is estimated 
to be small, in contrast to peripheral collisions. The data is not corrected for this effect, introducing 
a correction would reduce the scaled variance by less than 5%. 

In the NA49 experiment the number of target participants N p arg is not fixed, it fluctuates sub- 
stantially (see figure 181). These fluctuations might contribute even to multiplicity fluctuations in the 



forward hemisphere and increase the scaled variance [14, 16] 
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Figure 6: Left: multiplicity distribution of negative hadrons in central (C < 1 %) Pb 4- Pb collisions from 20A 
(top) to 158A GeV (bottom). The red line indicates a Poisson distribution with the same mean multiplicity. 
Right: ratio of the measured multiplicity distribution over Poisson distribution. It is visible that the measured 
distribution is significantly narrower than the Poissonian one for all energies. Note that the multiplicity 
fluctuations are not corrected for the finite centrality bin width. 
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Figure 7: Centrality dependence of co(h + ) (top), (o(hr) (middle) and co(h ± ) (bottom) for Pb+Pb collisions 
at different energies. C < 1% corresponds to the most central collisions. The shown centrality range of 
C < 10% corresponds approximately to number of projectile participants Np WJ > 160. Only statistical errors 
are shown. g 
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Figure 8: Fluctuations in the number of target participants for a fixed number of projectile participants in 
the UrQMD and HSD [g] models. 



Model calculations (HSD, UrQMD) [ 15] show that the scaled variance of target participant fluctu- 
ations co(N ! p rs ) = V ^^r g J is small for central collisions, and thus the corresponding increase of 
the scaled variance of multiplicity distributions is also expected to be small. 

In order to minimize the contribution due to both Veto calorimeter resolution and Np rg fluc- 
tuations further analysis was done for very central collisions (C < 1%, which corresponds to 
N p P roi > 193). 

The systematic error is estimated taking into account the influence of the track impact pa- 
rameter cut (±2%), the cut on electron rejection (±2%), the effect of finite resolution of the Veto 
calorimeter (—5%) and possible remaining participant fluctuations due to the finite width of the 
selected centrality bin (—5%). The total systematic error is estimated to be o sys %. 

Model calculations (Venus and UrQMD) show that for higher energies (80A and 158A GeV) 
the details of centrality selection have a significant effect on multiplicity fluctuations. When the 
centrality is determined by fixing the number of projectile spectators the scaled variance obtained 
by these models is up to 30% larger than the one obtained by fixing the energy in the Veto Calorime- 
ter using a detailed simulation of its acceptance. For lower energies (20A- 40A GeV) no significant 
difference of the scaled variance for different centrality selections was observed. 



4. Results and Discussion 

The results on energy dependence of multiplicity fluctuations in Pb + Pb collisions are shown 
in figure |9| It can be seen that at all energies the scaled variance for positively and negatively 
charged hadrons is smaller than 1, the value for a Poissonian distribution. 

A model independent comparison of results at different energies is not possible because the 
acceptance is different for each energy (figure ||). 
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Figure 9: Energy dependence of multiplicity fluctuations in the most central (C < 1%, Venus calculations: 
C < 2%) Pb + Pb collisions in comparison to string hadronic model (Venus and UrQMD) predictions for 
Co(h + ) (top left), co(h~) (top right) and o(/i ± ) (bottom). Horizontal bars indicate systematic errors of data 
points. 



10 



Energy Dependence of Multiplicity Fluctuations 



Benjamin Lungwitz 



The scaled variance for all charged hadrons is larger than the scaled variance for positively 
and negatively charged hadrons. This is probably due to the effect of resonance decays which 
correlated final state hadrons and thus leads to an increase of fluctuations. This effect is largest 
for all charged hadrons as the majority of resonances decay into one positively and one negatively 
charged daughter. It is smallest for negatively charged hadrons as only a very few resonances decay 
into two negatively charged hadrons (e.g. A — > %~ + p~). 

The multiplicity fluctuations of all charged hadrons in Pb+Pb collisions at the top SPS energy 
have been studied previously by WA98 A direct comparison of scaled variance of multiplicity 
distributions obtained in NA49 to the results of WA98 is not possible because of different exper- 
imental acceptance and centrality selection. The WA98 results were obtained in a pseudorapidity 
region around midrapidity (2.35 < f] < 3.75), whereas the NA49 results are obtained in the for- 
ward hemisphere (1.08 < y(n) < 2.57). The centrality of the collisions in WA98 was determined 
by fixing transverse energy in a calorimeter located in the pseudorapidity interval 3.5 < r\ < 5.5. 
In contrast in NA49 centrality is fixed using the energy in the forward Veto calorimeter (see sec- 



tion 2. 1 ) 



4.1 String- Hadronic Models 



For a comparison of data to string hadronic models (Venus 4.12 [ 17 ] and UrQMD 1.3 [18]) the 
model predictions were calculated in the data acceptance. The centrality selection in the models 
is done in a similar way as in the data by a detailed simulation of the acceptance of the Veto 
Calorimeter. 

The multiplicity distributions for negatively charged hadrons in data and models are compared 
in figure [K| In comparison to data the multiplicity distributions in Venus are shifted towards lower 
multiplicities, whereas in UrQMD they are shifted towards higher multiplicities. 

Figure |9| shows a comparison of the scaled variance between data and the predictions of the 
Venus and UrQMD model. The Venus model overpredicts scaled variance for central Pb + Pb 
collisions at most energies. In general the scaled variance obtained by the UrQMD model is in 
agreement with data. Only for a>(h + ) and (£>(hr) at 158A GeV as well as for ft)(/z ± ) at 80A and 
158A GeV the differences are significant. 

4.2 Acceptance Dependence 

Only a fraction of all produced particles is used for scaled variance determination in the experi- 
ment. A simple parametrisation of the acceptance dependence of ft) would be useful for comparison 
of different energies and a comparison to models. In Appendix |A| a derivation of a simple formula 
describing the dependence of the scaled variance on the acceptance is given. It is obtained under 
the assumption that all particles are produced independently in momentum space. If a fraction p of 
all particles N tot is accepted, the scaled variance ft) is related to the scaled variance (Ot ot in \% as: 



ft) = Wa J^J = l+p ((Otot — 1) (see equation |A.7| ). 



In figure 1 1 the fluctuations of negative hadrons for the five analysed energies are plotted as a 



function of the fraction of accepted tracks. If the simple acceptance dependence A.7 works and if 
the scaled variance in An is similar for all energies, the data points for all five energies should lie on 
a common straight line which starts at ft) = 1 for p = 0. It can be seen that the data seem to follow 
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Figure 10: Multiplicity distributions of negatively charged hadrons in the most central (C < 1%) Pb + Pb 
collisions from 20A to 158A GeV; data - black circles, Venus model - dashed red line and UrQMD - solid 
green line. Note that the multiplicity fluctuations are not corrected for the finite centrality bin width. 
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Figure 11: Acceptance dependence of scaled variance for negatively charged hadrons. The data points 
correspond from left to right to the energies of 20A to 158A GeV. Horizontal bars indicate systematic errors 
of data points. 



approximately the predicted scaling. A more detailed study of rapidity and transverse momentum 
dependence of multiplicity fluctuations is in progress, but first results indicate that the acceptance 
scaling seen in figure [Tl| should be considered only as a rough approximation. 
For co(h + ) this scaling does not work so well probably because a larger fraction of resonance de- 
cays (e.g. A ++ ) into two positively charged daughters, introducing more correlation in momentum 
space. The situation is even worse for co^) because both daughters of a decay into two unlike sign 
charged hadrons (e.g. p— > n + + 7i~) are used for the fluctuation analysis. Therefore acceptance 
scaling plots are shown only for (o(h~). 

4.3 Statistical Models 



Predictions of the statistical hadron-resonance gas model [19] on energy dependence of the 
scaled variance of negatively charged hadrons are shown in figure [12|. In this model multiplicity 
fluctuations in An were calculated for the limit of V — > °° using both canonical and grand-canonical 
ensembles. The effects of resonance decays and quantum statistics are taken into account. The 
calculated value ft) for the canonical ensemble is much lower than for the grand-canonical one in 
the SPS energy regime. In the canonical ensemble fluctuations are suppressed due to conservation 
laws of electric charge, baryon number and strangeness by about a factor of two. Preliminary 
calculations show a further suppression of fluctuations in the micro-canonical ensemble due to 
additional energy conservation. 



Figure 13 shows data on multiplicity fluctuations in comparison to the model calculations. 



The scaled variance in the limited acceptance was extrapolated using equation A..7 from the 



model results in full acceptance. As discussed above equation [A.7| may be considered as a reason 



able approximation only for co(h ). Use of equation A.7 introduces an uncertainty in the compar- 
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Figure 12: Energy dependence of the scaled variance of negatively charged hadrons in full acceptance in 
statistical hadron-resonance gas model [|l9[ for grand-canonical (GCE) and canonical (CE) ensemble. The 
"Final" lines show results after resonance decays. 
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Figure 13: Energy dependence of multiplicity fluctuations of negatively charged hadrons in Pb + Pb colli- 
sions in comparison to canonical and grand-canonical statistical hadron-resonance gas model [jig]. Horizon- 
tal bars indicate systematic errors of data points. 
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ison of data and model. Note also that in the model calculations the volume of the fireball is fixed, 
but in data the energy in the Veto calorimeter. 

The comparison in figure |13| shows that multiplicity fluctuations in data are much lower than 
predicted by the grand-canonical model. This is an indication of suppression of fluctuations in rel- 
ativistic gases due to conservation laws. The predictions of the canonical ensemble are also higher 
than the data. Energy- momentum conservation and the finite volume of hadrons are expected to 
cause an additional suppression of fluctuations. 



5. Summary 

Multiplicity fluctuations in the most central Pb + Pb collisions at 20A, 30A, 40A, SO A and 
158A GeV were studied by NA49 at CERN SPS. At all energies multiplicity distributions for pos- 
itively and negatively charged hadrons are significantly narrower than a corresponding Poissonian 
distribution. 

In comparison to data the multiplicity distributions in Venus are shifted towards lower mul- 
tiplicities, whereas in UrQMD they are shifted towards higher multiplicities. The Venus model 
overpredicts scaled variance for central Pb + Pb collisions at most energies, whereas the scaled 
variance obtained by the UrQMD model is in approximate agreement with data. 

Non-monotonic energy dependence of fluctuations, which might be related to the onset of 
deconfinement or a critical point, is not observed in data. 



The hadron-resonance gas model [19] in grand-canonical formulation overpredicts the mea- 
sured scaled variance. Introduction of the material conservation laws (canonical model) improve 
agreement with the data. 

The presented results may be the first experimental observation of the effect of conservation 
laws in relativistic gases. 
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A. Acceptance Dependence of Scaled Variance 

Let us assume that particles are emitted independently in momentum space. 
For a given total number of produced particles N tot , the probability to measure a number of particles 
N in a fixed acceptance p is a binomial distribution. 

P( N \ Ntot ) =P Ntot ( N ) = (A.l) 



15 



Energy Dependence of Multiplicity Fluctuations 



Benjamin Lungwitz 



Consequently the mean number of accepted particles is 

<N>=N tot p (A.2) 
and the variance of the binomial distributions is 

Var(N)=N totP (\-p) (A3) 

If one allows the total number of produced particles to vary according to distribution P(N tot ), the 
probability to measure a number of particles ./V is 

P{N) = J?P{N\N tol )P(N tot ) (A.4) 

N,„, 

The mean number of particles the in acceptance is 

<N>=p<N m > (A.5) 

The variance of N is given by: 

Var{N) =< Var{N\N tot ) > +Var{< N\N tot >) 

=< Var{N\N tot ) > +Var{pN tot ) (A.6) 
=< N tot > p{\ - p) +p 2 Var{N tot ) 

Finally the scaled variance is 

(0 = ^l = \+p{(O tot -\) (A.7) 

or 

co tot = ^P^ = l + (co-l)/ P (A.8) 
where co tot is the scaled variance of N tot . 
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